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Decomposition Products of Phosphine Under Pressure: PH2 Stable and Superconducting? 
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Evolutionary algorithms (EA) coupled with Density Functional Theory (DFT) calculations have been used to 
predict the most stable hydrides of phosphorous (PH„, n — 1 — 6 ) at 100, 150 and 200 GPa. At these pressures 
phosphine is unstable with respect to decomposition into the elemental phases, as well as PH 2 and H 2 . Three 
metallic PH 2 phases were found to be dynamically stable and superconducting between 100-200 GPa. One of 
these contains five formula units in the primitive cell and has C2/m symmetry (5FU-C2/m). It is comprised of 
ID periodic PH 3 -PH-PH 2 -PH-PH 3 oligomers. Two structurally related phases consisting of phosphorous atoms 
that are octahedrally coordinated by four phosphorous atoms in the equatorial positions and two hydrogen 
atoms in the axial positions (/4/mmm and 2FU-C'2/m) were the most stable phases between ~160-200 GPa. 
Their superconducting critical temperatures (Tc) were computed as being 70 and 76 K, respectively, via the 
Alien-Dynes modified McMillan formula and using a value of 0.1 for the Coulomb pseudopotential, p.*. Our 
results suggest that the superconductivity recently observed by Drozdov, Eremets and Troyan when phosphine 
was subject to pressures of 207 GPa in a diamond anvil cell may result from these, and other, decomposition 
products of phosphine. 

PACS numbers: 


INTRODUCTION 

The pressure variable can be used to synthesize materi¬ 
als with unique stoichiometries, and properties that would 
not be accessible otherwise IBil resulting in chemically 
different bonding motifs that are not observed at ambient 
conditions II2I EJ. Recent experimental work by Drozdov, 
Eremets and Troyan has shown that hydrides of the p-element 
sulfur are superconducting with a critical temperature, 
of 203 K at 150 GPa jO). This work was inspired by first- 
principles calculations suggesting that at P > 130 GPa H2S 
would transform to a phase that has a Tc of 80 K above 
160 GPa fTOl. The much higher critical temperature observed 
experimentally led to the suggestion that under pressure hy¬ 
drogen sulfide decomposes into a hydride with a stoichiome¬ 
try that is not stable at 1 atm, and this new phase is respon¬ 
sible for the remarkably high Tc M- Indeed, first principles- 
calculations im confirmed that the Tc (computed using the 
Alien-Dynes modified McMillan equation) of a novel Im3m- 
H3S phase matched exceedingly well with what was observed 
by Drozdov and co-workers. A plethora of theoretical calcu¬ 
lations have verified that the H3S stoichiometry is indeed the 
most stable at these pressures, analyzed the factors contribut¬ 
ing to Tc, examined anharmonicity as well as the isotope ef¬ 
fect, and computed Tc for the related compounds selenium hy¬ 
dride and tellurium hydride lfT2l4T^ . These studies highlight 
that new hydrogen-rich materials can be attained under pres¬ 
sure, some of which may have a Tc higher than those previ¬ 
ously thought possible for Bardeen-Cooper-Schrieffer (BCS)- 
type superconductors ll20l . and have invigorated the search 
for high-temperature superconductivity in high-pressure hy¬ 
drides. Eor example, recently the stability of hydrides with 
novel stoichiometries containing a group 15 element such as P, 
As, and Sb have been studied computationally, and the super¬ 
conducting properties of select systems were examined. El. 

Recently, new exciting experiments have revealed that 


pressure-induced high-temperature superconductivity may be 
found in other hydrides of a p-block element. Drozdov, 
Eremets and Troyan measured superconductivity in phos¬ 
phine, PH3, which was liquefied in a diamond anvil cell and 
subsequently compressed El. Resistance measurements re¬ 
vealed a Tc of 30 K and 103 K at 83 GPa and 207 GPa, respec¬ 
tively. Structural information on the superconducting phases 
was not provided. The possible existence of a high tempera¬ 
ture superconductor containing phosphorus and hydrogen mo¬ 
tivated us to examine the structural landscape of these ele¬ 
ments combined under pressure using the evolutionary algo¬ 
rithm XtalOpt. Similar to what was found for H„S under 
pressure, it is likely that the observed superconducting prop¬ 
erties do not arise from the hydride with the “classic” 1 atm 
stoichiometry, but instead decomposition products of phos¬ 
phine, such as PH2. Between the pressures of 100-200 GPa 
phosphine is thermodynamically less stable than PH2 and H2, 
and several superconducting PH2 phases are metastable with 
respect to solid hydrogen and phosphorous. The Tc computed 
for the PH2 structures via the Alien-Dynes modified McMil¬ 
lan equation are significantly larger than those expected for 
pure phosphorus, approaching the experimental values mea¬ 
sured in “phosphine”. Our findings therefore suggest that PH2 
may be another hydrogen-rich BCS-type superconductor. 

COMPUTATIONAL DETAILS 

A priori crystal structure prediction calculations were car¬ 
ried out using the open-source evolutionary algorithm (EA) 
XtalOpt Release 8 and 9 ll^ 1^ that has previously been 
used to predict the structures of a variety of binary hydrogen- 
rich phases under pressure ll25H29l . EA runs were carried out 
on the PH3 stoichiometry at 100, 150, and 200 GPa employing 
simulation cells with 1-6 formula units (EU) at 100 GPa and 2- 
3 EU at 150 and 200 GPa. In addition, structure searches were 
performed on the PH„, n = (1,2,4-6), systems at 100, 150, and 
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200 GPa using cells with 2-3 formula units, unless otherwise 
noted in Tables S1-S6 in the Supplementary Information (SI). 
Duplicate structures were detected via the XtalCOMP EqI 
algorithm. The lowest enthalpy structures from each search 
were relaxed in a pressure range from 100-200 GPa. 

Geometry optimizations and electronic structure calcula¬ 
tions were performed by using density functional theory 
as implemented in the Vienna Ab-Initio Simulation Pack¬ 
age (VASP) versions 5.2 and 5.4.1 OTl . with the gradient- 
corrected exchange and correlation functional of Perdew- 
Burke-Ernzerhof (PBE) ll^ . The projector augmented wave 
(PAW) method ll^ was used to treat the core states, and a 
plane-wave basis set with an energy cutoff of 700 eV was em¬ 
ployed. The His and P 3s/3p electrons were treated explic¬ 
itly in all of the calculations, using the POTCARs for PAW- 
PBE H and PAW-PBE P available in the potpaw-PBE.52.tar.gz 
file from the VASP repository. The fc-point grids were gener¬ 
ated using the E-centered Monkhorst-Pack scheme, and the 
number of divisions along each reciprocal lattice vector was 
chosen such that the product of this number with the real lat¬ 
tice constant was 30 A in the structure searches, and 40-50 A 
otherwise. The A:-meshes and energy cutoffs used resulted in 
enthalpies that were converged to within 1 meV/atom. It is 
likely important to employ density functionals approximating 
the effects of van der Waals (vdW) interactions for molecu¬ 
lar solids containing p-block elements near ambient pressure. 
However, it has been shown that at higher pressures the ef¬ 
fect of vdW interactions for these types of systems becomes 
negligible. Eor example, by ~40 GPa dispersion forces were 
expected to have a minimal effect on the structural parameters, 
and as a result the properties, of CO 2 ll34l . Eor this reason we 
expect the effect of vdW interactions on the results presented 
herein, which were obtained between 100 and 200 GPa, to 
be small. The Nudged Elastic Band (NEB) method llTSl was 
used to construct a reaction pathway between the lA/mmm 
and C2 jm PH 2 phases at 200 GPa. 

Phonon calculations were performed using the Quantum 
Espresso (QE) 13^ program to obtain the dynamical matrix 
and electron-phonon coupling (EPC) parameters. In the QE 
calculations, the H and P pseudopotentials, obtained from the 
QE pseudopotential library, were generated by the method 
of Trouiller-Martins with Is^ and 3s^3p^ valence configu¬ 
rations, respectively, along with the PBE generalized gradi¬ 
ent approximation. Plane-wave basis set cutoff energies were 
set to 80 Ry for all systems. The Brillouin-zone sampling 
scheme of MethfesseTPaxton using a smearing of 0.02 Ry and 
22 X 22 X 22 k-point grids were used for all 100 GPa calcu¬ 
lations of the PH 2 2PU-C'2/m and lA/mmm structures. At 
150 and 200 GPa we used k-point grids of 24 x 24 x 24 and 
30 X 30 X 30 for the PH 2 2PU-C'2 jm and lA/mmm structures, 
respectively. Eor the PH 2 5EU-C'2/m structure a 12 x 12 x 12 
k-point grid was used. Density functional perturbation theory 
as implemented in QE was employed for the phonon calcula¬ 
tions. The EPC matrix elements were calculated using 2x2x2 
g'-meshes for all of the structures at 100 GPa, as well as for 
PH 2 5PU-C'2/m at 150 GPa, and the phosphorus structures at 


all pressures. At 150 and 200 GPa 4x4x4 and 5x5x5 
g'-meshes were used for the PH 2 2PU-C'2/m and lA/mmm 
structures, respectively. The EPC parameter. A, was calcu¬ 
lated using a set of Gaussian broadenings in steps of 0.005 Ry 
from 0-0.300 Ry. The broadening for which A was converged 
to within 0.05 was between 0.015 and 0.040 Ry for all struc¬ 
tures. The critical superconducting temperature, Tc, has been 
estimated using the Alien-Dynes modified McMillan equation 
Jill as. 


T-T ^log 

T.= —exp 


1.04(1-f A) 


A-Ai*(l + 0.62A) 


( 1 ) 


where wiog is the logarithmic average frequency and is 
the renormalized Coulomb potential, often assumed to be be¬ 
tween 0.1 and 0 . 2 . 


RESULTS AND DISCUSSION 


PH„: Pressure Induced Decomposition of PH3 

Because the phases that PH 3 adopts under pressure are un¬ 
known, we have carried out evolutionary searches to predict 
the global minima configurations up to 200 GPa. A detailed 
analysis will be published elsewhere, but the coordinates for 
the most stable structures found at P > 100 GPa are pro¬ 
vided in the SI. Briefly, the 100 GPa phase is made up of 
layers of PH and hydrogen, whereas the 150 and 200 GPa 
structures consist of 1-dimensional (ID) PH 3 -PH 3 networks 
that resemble those present in the most stable PH 2 struc¬ 
ture at 100 GPa discussed below. The computed enthalpies 
of formation, AHf, of PH 3 with respect to the most stable 
structures of solid hydrogen and phosphorus ll40H44ll 

showed that the classic 1 atm stoichiometry is not thermody¬ 
namically stable at the pressures considered herein, and po¬ 
tentially even at lower pressures. The AHf for the reaction 
|H 2 (s) H- P(s) —>■ PH 3 were computed to be 58.0, 45.2, and 
42.6 meV/atom at 100, 150, and 200 GPa, respectively. 

Nonetheless, it is convenient to plot the AP^xn of the hy¬ 
drides of phosphorous, PH„ with n = 1 — 6, identified via 
our evolutionary searches for the reaction PH„ —(^^^)P -f 
(f )PH3 (where n = 0 - 3) and PH„ ^ {^)H 2 + PH3 
(where n = 3 — 6 ), as shown in the left and right hand sides 
of Pig. [TJa), respectively. Eor each of these plots a convex 
hull can be drawn, and the structures whose APrxn lie on the 
hull correspond to the most stable decomposition products of 
phosphine, or the most stable products arising from the reac¬ 
tion of phosphine with hydrogen. At 100, 150 and 200 GPa 
only PH2 is thermodynamically preferred over P/PH3, and 
the reaction PH3 —>PH2 + |H2 is exothermic. At 100 and 
150 GPa the reaction PH3-I-H2 —>PH5 is exothermic. Because 
PH5 becomes unstable by 200 GPa we focused our analysis on 
PH2 (the coordinates for PH5 are provided in the SI). 
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FIG. 1: (a) AiFrxn for the reaction: (left) PH„ ^ (5^)P+(f)PH3 
(where n = 0 — 3) and (right) PH„ —>• (^^)H 2 + PFIa (where 
n = 3 — 6) at 100, 150, and 200 GPa. The convex hulls are given by 
the dashed lines, (b) The enthalpies as a function of pressure for the 
5FU-C'2/m and 74/mmm phases of PH 2 relative to 2FU-C'2/m 
between 100-200 GPa. The dashed line indicates the transition pres¬ 
sure below which the 5FU-C'2/m phase has the lowest non-ZPE cor¬ 
rected enthalpy. The open symbols provide the relative enthalpies in¬ 
cluding the ZPE correction. At 150 GPa the ZPE-corrected enthalpy 
of 5PU-C'2/m (not shown on the plot) is 45.8 meV/atom higher than 
that of the 2FU-C'2/m structure. 


Evolutionary searches at 100, 150 and 200 GPa identified 
three unique low-enthalpy PH 2 phases, whose relative en¬ 
thalpies are provided in Fig. [^b). A C2/m symmetry sys¬ 
tem whose primitive cell contained five formula units, 5FU- 
C2lm, had the lowest enthalpy excluding zero point energy 
(ZPE) contributions below 153 GPa. A two formula unit cell 
with C2/m symmetry (2FU-C'2/m) and an lA/mmm struc¬ 
ture became preferred at higher pressures. Even though these 
species are thermodynamically unstable relative to solid phos¬ 
phorous and H 2 (by at least 38.7, 29.4, and 31.9 meV/atom 
at 100, 150 and 200 GPa), they are all dynamically stable at 
100, 150 and 200 GPa. Moreover, the formation of PH„ with 
n = 4 — 6 from PH 2 and H 2 at 100, 150 and 200 GPa, and the 
decomposition of PH 2 into H 2 and PH at 100 and 150 GPa is 
computed as being endothermic. 

The coordination and arrangement of the phosphorus atoms 
in 5FU-C'2 /to is quite different than in the 2F\J-C2/m and 
lA/mmm phases. The former consists of PH 3 -PH-PH 2 -PH- 
PH 3 oligomers extending along the c-axis, which are ID pe¬ 
riodic along the a-axis, J^[PH 3 -PH-PH 2 -PH-PH 3 ], see Fig. 
|^a,b). The PH 3 -PH-PH 2 -PH-PH 3 unit possesses an inversion 


center of symmetry, with the H 3 P-PH and HP-PH 2 distances 
measuring 2.191 and 2.167 A, respectively, at 100 GPa. The 
phosphorous atoms comprising the PH and PH 2 units are co¬ 
ordinated to four other phosphorous atoms forming a square 
net. Such structural motifs are not uncommon for compressed 
phosphorous, particularly in the pressure range between 100 - 
200 GPa. For example, a number of the high pressure phases 
of elemental phosphorus, such as the A17 and A7 phases, have 
been described as Peierls distortions of a simple cubic lattice 
Hm |45l |46l. Beyond the A7 phase of black phosphorus, a 
simple cubic lattice is adopted at ~ 11.1 GPa that is stable up 
to 137 GPa Enumiiii. After a relatively brief transition 
through a suspected incommensurate phase BtI l48l . simple 
cubic phosphorus transitions into a simple hexagonal phase 
II 42 II followed by a body centered cubic phase adopted at pres¬ 
sures above 262 GPa ||43]| . 



FIG. 2: (a) A 3x3x2 supercell of the 100 GPa 5FU-C2/m struc¬ 
ture. Phosphorous atoms are orange, hydrogen atoms are white, (b) 
A fragment of the ^[PH 3 -PH-PH 2 -PH-PH 3 ] motifs that comprise 
5FU-C'2/m. The values provided in black are P-P distances (A) 
and the red values are bond angles at 100 GPa. (c) Electronic band 
structure along with the total and site projected electronic DOS of 
5FU-C'2/m at 100 GPa. Because the PDOS for the H and P atoms 
were relatively independent upon the local atomic environment, we 
have not decomposed the PDOS into contributions from distinct H 
and P atoms. Ff is set to zero. 

The phosphorous atoms in the -PH 3 units at the end of 
the oligomer are nearly perfectly octahedrally coordinated to 
three hydrogen and three phosphorous atoms, with the H-P-H 
angles measuring 87.4 and 92.5° and the P-P-P angles measur¬ 
ing 90° at 100 GPa. Similarly, the phosphorous atoms com¬ 
prising the -PH 2 - units at the center of the motif also assume 
octahedral coordination with four phosphorous atoms on the 
equatorial and two hydrogen atoms on the axial positions. The 
phosphorous atoms within the -PH- segment of the oligomer, 
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on the other hand, are nearly trigonal bipyramidal. The ax¬ 
ial positions are occupied by two phosphorous atoms, and the 
equatorial positions are filled by two phosphorous and one hy¬ 
drogen atom (ZPPP = 131.1°, ZPPH = 113.7 and 115.2°). 
At this pressure all of the P-H distances fall between 1.430- 
1.444 A, which is close to the typical experimental P-H bond 
length of 1.437 A ll49l at 1 atm. The P-P distances fall in 
the range of 2.167-2.191 A within the oligomer, and measure 
2.151 A between oligomers, which is comparable to the dis¬ 
tance we calculate for elemental phosphorous at this pressure, 
2.149 A. 

Fig.|^c) shows that at 100 GPa 5 FU-C' 2 /m-PH 2 is metal¬ 
lic. The occupied density of states (DOS) is nearly free elec¬ 
tron like. But because the Fermi level, Ep, falls near a pseu¬ 
dogap the DOS at Ep, g{Ep) = 0.045 eV“^/valence electron, 
is lower than the value computed for a free electron gas with 
the same bandwidth, 0.063 eV“ ^/valence electron. The metal- 
licity is due primarily to hydrogen states with s-character and 
also phosphorous states with p-character. The occupied hy¬ 
drogen and phosphorous states hybridize, indicative of cova¬ 
lent bonding. The electron localization function (ELF), pro¬ 
vided in the SI, shows regions of high ELF along the P-P 
and P-H bonds, and no lone pairs are observed, as would 
be expected based upon valence-shell electron-pair repulsion 
(VSEPR) theory, which assumes that the geometry a molecule 
assumes minimizes the repulsion between valence shell elec¬ 
tron pairs. A Bader analysis revealed electron transfer from 
phosphorous to hydrogen, with the average charge on the hy¬ 
drogen atoms being -0.34 at 100 GPa. The charges on the 
phosphorous atoms in the -PH 3 , -PH 2 -, and -PH- units were 
computed as being H-L12, h- 0.70 and h- 0.25, respectively. At 
150 GPa the average Bader charges become -0.37 (H), H-L15 
(P in -PH 3 ), h- 0.75 (P in -PH 2 -) and h- 0.30 (P in -PH-). Phonon 
calculations showed that this phase is dynamically stable at 
150 GPa. 

Above ^160 GPa two structurally related PH 2 phases, 
2PU-C'2/m (Pig. 1^ and lA/mmm (Pig. |^, become more 
stable than 5PU-C2 /to. The differences in the enthalpies 
of these two compounds at 150 and 200 GPa are less than 
1 meV/atom, but in between these pressures the non-ZPE 
corrected enthalpy of the lA/mmm phase is slightly lower 
than that of 2PU-C'2/m. The existence of nearly isoenthalpic 
hydrogen-rich phases under pressure in the DPT calculations 
is not unprecedented. We have previously found nearly isoen¬ 
thalpic but distinct structures in our theoretical studies of hy¬ 
drides containing an alkali metal or alkaline earth metal un¬ 
der pressure. This includes two LisH phases that were com¬ 
puted as being more stable than Li and LiH above 50 GPa lIZTl . 
five unique CSH 3 structures that were preferred over CsH and 
H 2 between 30-150 GPa ESl and three BaHg phases that had 
nearly the same enthalpies around 100 GPa ED- 

In both the C2/m and lA/mmm phases of PH 2 each 
phosphorous is octahedrally coordinated by four phosphorous 
atoms in the equatorial positions, and two hydrogen atoms in 
the axial positions. The difference between the two struc¬ 
tures is that whereas in lA/mmm all of the atoms assume 




FIG. 3: (a) (left) A conventional unit cell and (right) a 1x2x2 su¬ 
percell of the 200 GPa 2FU-C'2/m phase. Phosphorous atoms are 
orange, hydrogen atoms are white, (left) The canting angle with re¬ 
spect to a line perpendicular to the plane containing the phosphorous 
atoms (10.1°) is shown by the arrows, (right) The octahedra around 
the central plane of phosphorous atoms are highlighted by teal and 
gray, (b) Electronic band structure along with the total and site pro¬ 
jected electronic DOS of 2FU-C'2/m at 200 GPa. Ep is set to zero. 


the ideal octahedral angles, in 2FU-(72/m the hydrogens are 
canted with respect to a line that lies normal to the phospho¬ 
rous square net. At 100 GPa the canting angle is negligi¬ 
ble (<0.1°), and the P-P and P-H bond lengths in the two 
structures are nearly identical (Table|^. At 160 and 200 GPa, 
however, the canting angle increases to 9.8 and 10.1°, re¬ 
spectively. Even though the nearest neighbor P-P/P-H dis¬ 
tances in lA/mmm are slightly smaller than in 2FU-C'2/m 
at 160 and 200 GPa, the volume of the former is somewhat 
larger than that of the latter due to the tilting of the octahe¬ 
dra. The more efficient packing of the octahedra between the 
layers, which is a consequence of the canting, can most easily 
be seen by comparing the P-P distances between the phos¬ 
phorous atoms comprising different layers. At 100 GPa the 
separation is identical, 2.983 A, resulting in comparable vol¬ 
umes. By 160 GPa, however, this distance becomes 2.757 and 
2.695 A for lA/mmm and 2PU-C'2/m, respectively, giving 
rise to the sudden difference in volume of 1 . 2 % between the 
two phases. The PV contribution to the enthalpy is smaller in 
the more compact 2PU-C'2 /to phase, but the electronic con¬ 
tribution more strongly favors the more symmetric lA/mmra 
structure (see the SI), such that lA/mmm is preferred by a 
few meV/atom for P ~155-195 GPa. By 200 GPa, however, 
the two phases become isoenthalpic. 

Similar to 5PU-C2 /to, regions of high ELE in 2EU-C'2/m 
and lA/ramm are found only along the P-P and P-H bonds 
and the Bader charges in Tablej^are indicative of charge trans¬ 
fer from phosphorous to hydrogen. Both of these phases are 
good metals with g{Ep)= 0.045-0.048 eV“^/valence electron 
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TABLE I: Phosphorus-phosphorus, phosphorus-hydrogen and hydrogen-hydrogen distances for the 2FU-C'2/m and I'i/mmm PH 2 phases. 
P-H(l)/H-H(l) are the first nearest neighbor, and P-H(2) are the second nearest neighbor distances. The average Bader charges per atom type 
are provided, along with the volume per unit cell. 


System 

Pressure P-H(l) P-H(2) 
(GPa) (A) (A) 

P-P 

(A) 

H-H(l) 

(A) 

H 

(e) 

P 

(e) 

Volume 

(A^/FU) 

2FU-G2/m 

100 

1.449 

2.165 

2.159 

1.527 

-0.26 

0.53 

13.885 


150 

1.436 

2.013 

2.111 

1.497 

-0.29 

0.57 

12.398 


160 

1.445 

1.760 

2.127 

1.510 

-0.27 

0.55 

12.025 


200 

1.439 

1.693 

2.093 

1.500 

-0.26 

0.52 

11.288 

lilmmm 

100 

1.449 

2.163 

2.157 

1.528 

-0.30 

0.60 

13.882 


150 

1.437 

2.012 

2.109 

1.494 

-0.28 

0.56 

12.401 


160 

1.435 

1.989 

2.101 

1.490 

-0.28 

0.56 

12.170 


200 

1.434 

1.896 

2.077 

1.487 

-0.26 

0.52 

11.359 



FIG. 4: Same as Fig. except for PH 2 at 200 GPa. 

/PPP=90°, ZHPP=90° and ZHPH=180°. In the figure on the top 
right a 2x2x 1 supercell is shown. 


at 200 GPa, which is close to the value computed for a free 
electron gas with the same bandwidth, 0.054 eV“^/valence 
electron. Their DOS plots display many features in common 
with the DOS calculated for 5FU-C'2/m: they are indica¬ 
tive of H/P hybridization and the character at the Fermi level 
is primarily hydrogen s-like, with substantial contributions 
from phosphorous p-states. Phonon calculations revealed that 
these two structures are dynamically stable at 100, 150, and 
200 GPa. 

Because of the structural similarity between the 2FU- 
C2/m and lA/mmm phases, we wondered if they may be 
connected via a low energy pathway. The NEB method was 
employed to find the barrier and transition state between the 
two structures. The computed barrier was ^270 meV/atom 
at 200 GPa, implying that these phases likely lie within two 
distinct wells on the potential energy surface. The transition 
state (whose coordinates are given in the SI) can best be de¬ 


scribed as the average of the two PH 2 phases. Within the NEB 
method, the positions of the individual atoms within the start¬ 
ing and ending structures must match as closely as possible 
to avoid large geometric changes that may result in unrealistic 
barriers. In order to do this the lA/mmm structure has to be 
reduced down to its primitive triclinic cell and then a 2 EU 
supercell in the crystallographic 6-direction direction must be 
constructed. In the transition state structure the perfect 90° 
square phosphorus nets described above have been distorted 
to diamonds with internal angles of 67 and 113° resulting in 
nearest and second nearest neighbor P-P distances of 1.885 
and 2.581 A (c.f. Table 1), and the distance separating the 
P-P layers is 2.941 A, which can be compared to 3.017 and 
2.963 A for the I4/minm and C2/m phases, respectively. 
Each phosphorus is [2 h- 4] coordinate to hydrogen. Similar to 
the two end point structures, the shortest two P-H distances 
fall between 1.45 and 1.50 A, which is slightly longer than 
the distance of a P-H bond in phosphine at ambient conditions, 
1.42 A II 50 I . The remaining four P-H distances fall within 1.7- 
1.8 A, and are different for each phosphorous atom in the PI 
symmetry 2PU transition state. 

We also wondered if tilting the octahedra to different angles 
would yield other structures not found in our evolutionary runs 
whose enthalpies were similar to lA/mmm and C2/m. In 
order to test this hypothesis we relaxed structures constructed 
to have canting angles of 45°, as well as geometries where 
hydrogen atoms were placed at random positions within the 
phosphorous lattice at 200 GPa. No unique low enthalpy 
structures were found. In addition, a closer inspection of the 
evolutionary runs showed that structures similar to those we 
constructed by hand had been generated and optimized during 
the EA searches, but their enthalpies were significantly higher 
than 5EU-C'2/m, 2FU-C'2/m and M/mmm PH 2 . 

Because the ZPE may be large for extended systems con¬ 
taining light elements, we computed the enthalpic differences 
between these three PH 2 phases including the ZPE correc¬ 
tions at pressures where they were dynamically stable. The 
open symbols in Eig. [^c) show that the inclusion of the ZPE 
makes the lA/mmm and 2EU-C'2/m structures significantly 
more stable than 5EU-C'2/m at 150 GPa, with the C2lm be- 
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ing preferred by 6 meV/atom. Because the 2FU-C'2/m and 
lA/mmm phases are structurally similar, so is their phonon 
DOS, and their ZPE at 150 GPa is nearly identical. The 
phonon DOS of the 5FU-C'2/m phase, on the other hand, 
shows that the frequencies associated with the hydrogen vi¬ 
brations above 2000 cm“^ are larger in comparison to these 
two phases, resulting in a significantly larger ZPE at 150 GPa, 
which leads to the destabilization of this phase when com¬ 
pared to the other two structures. At 200 GPa the ZPE cor¬ 
rections favor the lA/mmm phase such that its zero-point 
corrected enthalpy is 25 meV/atom lower than that of 2EU- 
C2/m. The reason for this is that the C^jm structure has 
more high frequency modes above 1500 cm“^ (see Fig. |^. 
Thus, our results present another example of how the ZPE of 
phases with light elements can affect their relative enthalpies. 
It is beyond the scope of this work to consider anharmonic 
effects to the (free) energies (and Tc), but it may be that 
these also influence the order of stability of the phases stud¬ 
ied herein. To determine if these three PH 2 structures could 
potentially contribute to the superconductivity observed by 
Drozdov and co-workers in their compression of phosphine 
up to 207 GPa, their superconducting properties were investi¬ 
gated in further detail, as described in the following section. 

Superconducting Phases of PH2 

Above 5 GPa superconductivity has been observed in black 
phosphorous, the most stable allotrope at ambient conditions. 
The experimentally measured depends on the path taken 
in the pressure/temperature phase diagram. Recent experi¬ 
ments, and DFT calculations using the Alien-Dynes modified 
McMillan equation, both showed that for the simple cubic 
phase of phosphorous Tc decreases with increasing pressure 
15211531 . Computations revealed that the decrease in Tc above 
30 GPa could be explained by the increase in the phonon fre¬ 
quencies, and a /r* of 0.18 yielded a Tc of 5.8 K at 70 GPa 
II 53 I . which agrees well with the experimentally measured Tc 
of 4.3 K at 100 GPa IST). At higher pressures, up to 160 GPa, 
experiments suggested that Tc decreased below 4 K 15^ . In 
comparison, the computational methodology used herein cou¬ 
pled with a /r* of 0.18 resulted in a Tc of 6.4 K for simple 
cubic phosphorous at 100 GPa, and slightly higher tempera¬ 
tures for smaller values of the Coulomb pseudopotential, see 
Table |I^ We have also computed the parameters entering the 
Alien-Dynes modified McMillan equation for simple hexag¬ 
onal phosphorous at 150 and 200 GPa, so that they may be 
compared with those calculated for the aforementioned PH 2 
phases. As Table[n|shows, the average logarithmic frequency, 
wiog, of phosphorous increases with pressure and the electron- 
phonon coupling. A, decreases quenching the superconductiv¬ 
ity within the simple hexagonal phase. 

How would the incorporation of H atoms within the phos¬ 
phorus lattices, in the most stable PH 2 geometries, affect the 
uj\og, A and Tc? To answer this question we computed the 
superconducting properties of the previously discussed PH 2 


phases, at pressures where they were found to be dynamically 
stable, and the results are given in Table Fig. provides 
the phonon band structure and linewidths, Eliashberg spectral 
function, and plots the dependence of A on the frequencies 
for the 2FU-C'2 /to and lA/ramm phases at 200 GPa (analo¬ 
gous plots can be found in the SI for all of the other Tc val¬ 
ues computed herein). In pure compressed hydrogen it has 
been proposed that p* can range from 0.08 to 0.08911541 l55l. 
whereas values of 0.18 ||5T| have been employed for simple 
cubic phosphorus at 100 GPa. We therefore list critical tem¬ 
peratures computed for /r* of 0.1 and 0.18 in Table [H] How¬ 
ever, because the generally accepted value for the Coulomb 
pseudopotential is within 0.1 to 0.13, the text quotes Tc val¬ 
ues computed for a fj,* of 0.1. 

Because of the structural similarities between C2/m and 
lA/mmm at 150 GPa, their critical temperatures should be 
nearly identical, and Tc was calculated as being 51 K for the 
lA/mmm structure. The aforementioned canting of the octa- 
hedra results in geometric changes between these two struc¬ 
tures at higher pressures, which could potentially affect the 
superconducting properties. Indeed, at 200 GPa the electron 
phonon coupling of lA/mmm was calculated as being some¬ 
what larger than for C2lm, but wiog was nearly 175 K larger 
for C2/m because the vibrational modes arising primarily 
from hydrogen atoms are, on average, found at higher fre¬ 
quencies in this phase. The canting results in a larger number 
of hydrogen atoms that are found within a distance of 2 A of 
a given hydrogen in 672/771, thereby increasing the frequen¬ 
cies associated with their vibrations. For example, whereas 
lA/mmm has four H-H distances at 1.486 A, in C2/m the 
H-H measures are 1x1.497 A, 1x1.507 A, 2x1.585 A, and 
2 X 1.775 A. Because of the opposing effects of wiog and A on 
the critical temperature, the Tc values computed for the 672/m 
(76 K) and lA/mmm (70 K) phases at 200 K are very similar 
despite their structural differences. 

At 200 GPa the most significant contribution to A in the 
2FU-672/m and lA/mmm PH 2 systems stems from vibra¬ 
tional modes between 500-1750 cm“^ that involve hydrogen 
and phosphorus motions. The electron phonon coupling aris¬ 
ing from modes above 1750 cm“^, which consist primarily 
of the motion of hydrogen atoms with a small share from 
the phosphorous atoms, are also substantial components of 
the total A. Vibrations below 500 cm“^, which are the re¬ 
sult of phosphorus motions, play a small, but non-negligible 
role in the electron-phonon coupling. The computed Tc of 
the PH 2 phases are significantly larger than for monoatomic 
phosphorus because hydrogen increases both wiog and A. The 
electron phonon coupling is in-line with the values typically 
computed for compressed hydrogen-rich phases, 0.5-1.6 ll56l - 
l58l . but smaller than the values of 2.19 at 200 GPa and 2.69 
at 150 GPa computed for Im3m H 3 S ifTTl and CaHg 1591 . 
respectively. The phonon linewidths reveal that unlike CaHg 
I 59 II . where the electron phonon coupling was derived primar¬ 
ily from a single vibrational mode, a plethora of modes con¬ 
tribute towards A in PH 2 . 

The superconducting critical temperatures we compute for 
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TABLE II: Electron-phonon coupling parameter (A), logarithmic average of phonon frequencies (tuiog) and estimated superconducting critical 
temperature (Tc) for values of the Coulomb pseudopotential (fi*) of 0.1 and 0.18 for simple cubic (S.C.) phosphorous, simple hexagonal (S.H.) 
phosphorous, and the 5FU-C'2/m, 2EU-C2/m, and lA/mmm PH2 phases at various pressures. 


System 

Pressure (GPa) 

A 

(K) Tf 

= 0.1 

(K) (K) 

S.C.P 

100 

0.66 

521.7 

15.9 

6.4 

S.H. P 

150 

0.21 

575.6 

0.00 

0.00 

S.H. P 

200 

0.13 

671.8 

0.00 

0.00 

5FU-C'2/m 

100 

1.05 

655.1 

49.0 

32.1 

5FU-C'2/m 

150 

1.00 

798.1 

55.5 

35.2 

2FU-C'2/m 

200 

1.04 

1026.5 

75.6 

49.2 

I4/mmm 

150 

0.86 

946.2 

50.6 

28.3 

I4/mmm 

200 

1.13 

851.6 

70.4 

48.0 




FIG. 5: Phonon band structure, Eliashberg spectral function, 
o?F{uj), and the electron-phonon integral, A(cu), for 2FU-C'2/m 
and lAjmmm PH2 at 200 GPa. Circles indicate the phonon 
linewidth with a radius proportional to the strength. The highlighted 
sections of the a‘^F{uj) plots show the division into vibrational 
modes that are comprised of phosphorus vibrations (P-P, green), mo¬ 
tions of hydrogen and phosphorus atoms (P-H, pink), and mainly 
hydrogen vibrations (H-H, white). The contribution towards A for 
these divisions are 0.09 (9.3%), 0.50 (54.6%), and 0.33 (36.1%) for 
2FU-C'2/m and 0.13 (12.8%), 0.58 (54.9%), and 0.34 (32.3%) for 
M/mmm, for the green, pink, and white divisions, respectively. 

the lA/mmm and 2FU-C'2 /to phases of PH 2 are somewhat 
lower than the experimentally measured value of 103 K at 
207 GPa EB. However, it may be that the measured by 
Drozdov and co-workers results from a mixture of phases. 
For example, recent experimental work has suggested that 
the decomposition of LiH under pressure in a diamond anvil 
cell may lead to the formation of layers with different LiH„ 
stoichiometries II 6 OII . The experiments were inspired by the¬ 
oretical predictions ll^ . and observables computed via DFT 
calculations were employed to aid the interpretation of the ex¬ 
perimental results. Pepin and co-workers suggested that un¬ 
der pressure Li diffuses into the diamond anvil cell forming a 


LiHg layer at the diamond/sample interface, and an LiH 2 layer 
at the LiHg/LiH interface. Such mechanisms may also be im¬ 
portant for the pressure induced decomposition of phosphine, 
and it is only via comparison of the computed experimental 
observables for specific phases with the results obtained ex¬ 
perimentally that one can uncover which phases are formed 
under pressure, and the mechanisms underlying their forma¬ 
tion. A feedback loop between experiment and theory is inte¬ 
gral to advance our understanding of high pressure phenom¬ 
ena. 

The decrease in with pressure of all of the PH 2 phases 
considered herein correlates with the markedly lower critical 
temperature, 30 K, measured by Drozdov and co-workers at 
83 GPa. In future work we will focus on predicting the struc¬ 
tures of PH„ phases at pressures where vdW interactions may 
be important, and interrogate their bonding and superconduct¬ 
ing properties. 

CONCLUSIONS 

Recent experiments revealed that when PH 3 is compressed 
to 207 GPa, it becomes superconducting below 103 K ED. 
Our density functional theory (DFT) calculations have shown 
that at pressures of 100, 150 and 200 GPa PH 3 is ther¬ 
modynamically unstable with respect to decomposition into 
the elemental phases, as well as PH 2 and H 2 . Based upon 
the computed enthalpies other reactions that may occur un¬ 
der pressure are: PH 3 +H 2 —^PHg (100 and 150 GPa) and 
PH 2 —|H 2 H-PH (100 GPa). A priori crystal structure pre¬ 
diction has been used to identify three PH 2 phases that are 
dynamically stable in the pressure range of 100-200 GPa. 
All of these are calculated to be superconducting via the 
Alien-Dynes modified McMillan equation, suggesting that 
like H 3 S a, PH 2 may be another hydrogen-rich Bardeen- 
Cooper-Schrieffer (BCS)-type superconductor. 

Two PH 2 structures with C2/m and lA/mmm symmetry 
were computed to have a superconducting critical tempera¬ 
ture, Tc, of 76 K and 70 K, respectively, at 200 GPa. Between 
150-200 GPa the non-ZPE (zero point energy) corrected en¬ 
thalpies of these phases differed by only a few meV/atom, 
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and their structures both consisted of square nets of phospho¬ 
rous atoms. In the li/mmm structure the phosphorous atoms 
were octahedrally coordinated, with two hydrogen atoms in 
the axial positions. In the C'2/m phase the hydrogens are 
slightly canted from the ideal octahedral angles allowing the 
distance between the planes of phosphorous atoms to de¬ 
crease, and reducing the volume of the structure. A five for¬ 
mula unit phase with C^jm symmetry (5FU-C'2/m), which 
consists of 1-D PH 3 -PH-PH 2 -PH-PH 3 oligomers, was also 
identified in our evolutionary searches, and it has the lowest 
non-ZPE corrected enthalpy at 100 GPa. At this pressure it’s 
Tc is computed to be 49 K. Upon decreasing the pressure, the 
Tc of all of the PH 2 phases decreased. 

Our results provide another example of how pressure can 
lead to the formation of compounds with stoichiometries and 
properties that would not be predicted based upon our expe¬ 
rience at 1 atm. Comparison of the observables computed 
for these PH 2 phases with results obtained experimentally 
will unveil which phase, or mixture of phases, give rise to 
the superconducting properties observed by Drozdov and co¬ 
workers ED- 

Supplementary Information Available: Details of the 
structure searches, band structures, structural parameters, 
DOS plots, phonon DOS plots, absolute energies in Hartree, 
phonon convergence plots, and Nudged Elastic Band results. 
This material is available free of charge via the Internet at 
http://pubs.acs.org 
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